function curvDiff = compareCurvature( bwImage, pixelA, pixelB, kernelRadius )
%UNTITLED Summary of this function goes here
%   Returns the difference in curvature at two points as estimated by the 
%   Local Area Integral Invariant.  

       laii = getLaiiSignature(bwImage, kernelRadius);
       laiiA = laii(pixelA(1,1), pixelA(1,2) );
       laiiB = laii(pixelB(1,1), pixelB(1,2) );
       curvDiff = laiiA - laiiB;
end

